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CHAPTER I 


INTRODUCTION 


Introd u cj bipn; 

In recent years, much interest has been focussed on various 
problems that arise in cyclic staffing. Such problems involve 
the optimal scheduling of resources to meet demands where both 
resource availability and demand profile are cyclic. 

Many industries and public services operate around the 
clock, 7 days a week. This is the case in police and fire depart- 
ments, in telephone companies, in hospitals and in many heavy 
industries. Typical work schedules for employees of such services 
and industries normally consist of a cycle spreading over a few 
weeks and in which shift work and holiday periods alternate. 

These schedules consist of a cycle which repeats itself after a 
few weeks. Employees can therefore be seen as continuously 
going throiigh the same basic work-holiday pattern or rotating 
schedule. One rotating schedule is determined by the length of 
its cycle, by the length and frequency of the work and rest 
periods and by the way these alternate. 

The obcjective in cyclic staff scheduling problem is to 
minimize the total labor hours needed to meet demand requirements 
during different time periods. This problem is efficiently 
solvable when resources are continuously available during conse- 
cutive time periods. For an example of continuous availability, 
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let us consider one fundamental problem, called (5-7) cyclic 
staffing problem. In this problem, both resource availability 
and requirements are cyclic. The obj active is to minimize the 
number of workers needed to meet daily requirements, with the 
added provision that each worker be allowed two consecutive days 
off each week. 

Let X. denote the number of workers in shift j and b. 

O •>*“ 

denote the requirement during time period i. With thesen nota- 
tions the above problem can bo formulated asi 
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X 0 and integer for all o = 1> 7. 

i J 

Here the \;ork-shift is five consecutive days. If wo 
denote the 0-1 matrix by A, each column of A represents a possi- 
ble work shift and each row of A shows the shifts that are active 
on the corresponding day. That the resources arc continuously^ 
available is reflected in A, where each column contains exactly 
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one block of consecutigo ones (when viev/ed circulsrily) . This 
problem has been solved optimally by Bartholdi et al [6], 

In a more realistic model, resources are intermittently 
available. This may be due to, for example, limch and relief 
breaks for the workers. The following matrix is an hypothetical 
example showing intermittently available resources. 

100011 01 ! 

11000110 

01100011 

10110001 

110 11000 

01101100 

00110110 

00011011 

L 

This IS an (8 X 8) matrix and shift length is five. It 
can be observed that there is one break (when viewed circularity), 
i.e. two blocks of consecutive ones, in each row, 

Bartholdi [5] showed that the problem of cyclic staff ^ 
scheduling with intermittently available resources is a NP- 
complete problem. In view of this observation, the possibility 
of developing an exact algorithm for the above problem is remote. 
Heuristic approaches based on linear programming relaxation, 
network relaxation or lagrangean relaxation are available in 
literature. In all of them, some explicit use was made of the 
fact that A almost had the property of consecutive ones. 
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Bartholdi [ 5 ] developed an algorithm based on linear 
programming round-ofi by making no assumptions about the circu- 
larity of A, He gives a bound on the maximum possible deviation 
of his heuristic from the integer optimum solution. This bound 
is not good for the problems in which interniittent breaks are 
reasonably large in each siiift. More often than not, it is 
observed that the deviation of the solution obtained using 
Bartholdi's method from the integer optimum solution is almost 
equal to the bound given by him. 

In this dissertation, we consider the problem of scheduling 
with intermittently available resources and develop some methods 
for solving it. We develop three algorithms. One is LP-based 
heuristic algorithm. In this algorithm, a sequence of IP’s are 
solved and an improved version of the rounding-off procedure 
suggested by Bartholdi et al [6] is used after solving each LP 
to got a feasible solution to the cyclic staff scheduling problem. 
The other two algorithms are Branch and Bound algorithm and cutting 
plcine algorithm. These two are exact algorithms. Branch and Bound 
algorithm makes use ol sensitivity analysis to impose lower or 
upper bound constraints on any variable and LP-based heuristio 
is applied at every mnodo in the branching tree. The cutting 
plane algorithm makes use of the dual fractional cutting plane 
algorithm. Two types of cuts are introduced to reduce the total 
number of cuts needed. 
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All the suggested algorithms have been programmed and 
tested extensively on randomly generated cyclic staff schedul- 
ing problems. The results are quite encouraging and reported. 

The solutions obtained using LP-based heuristic algorithm are 
observed to be very close to the optimum solution for all the 
problems solved. The Branch and Bound algorithm is able to 
solve problems of size (96 x 96 ) optimally in a few minutes of 
CPU time. The performance of cutting plane algorithm was 
observed to be erratic. It could not solve problems of size 
more than (48 x 48). 

1 , 2 Brel imi nar ies;: 

Some notations and well-known concepts of cyclic staff 
scheduling problem are used throughout the thesis. For the 
sake of completeness they are given below, 

A 0-1 vector is said to be circular if its I's occur 
consecutively, where the first and last entries are considered 
to be consecutive. 

A matrix is called column circular (resp., row circular) 
if its columns (resp., rows) are circular. 

The term scheduling in this context, denotes the act of 
specifying how many resources should there be in different shifts 
so as to satisfy the requirements dinring different time periods 
with minimum cost. 
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Following notations are used in the thesis: 


m 

n 

A 

X 


a. 


10 ’ 


k 




c 

b 

0 

1 

L= 1 
r. I 


total number o£ time intervals to be scheduled, 
total nmber of different shifts or work schedules, 
an (n X n) matrix of O's and I’s. 
number of operators assigned to the Q-th work 
schedule . 

number of operators required for time interval i. 

^ 1 if time interval i is a work period in the j“th 
work schedule. 

0 otherwise . 

V 

cost of assigning a worker to shift 

maximum number of blocks of consecutive ones (when 

considered circularit^O m an^^ rov/ in matrix A, 

minimum number of ones in each column of A. 

assignment vector (n x l). 

cost vector (l x n). 

requirement vector (m x l). 

vector of zeros (m x l). 

vector of ones (l x n) . 

the greatest integer less than or equal to . 

the least integer greater than or equal to . 


1,3 Stat eme nt pf^ th-iL FJlSi’iPR* 

The problem of cyclic staff scheduling can be posed as 
follows s given resource requirements for each time period within 
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a planning horizon, find out a schedule of resources to meet 
these requirements at minimum cost. This can be formulated as 
an IIP as given below. 

n 

Min. I. c . X 
0=1 ^ ^ 

s.t . 

n 

I a. .x^ b , 1 = 1, . . . , m 

0=1 ^ ^ 

X 0, 3 = 1, . , . , n 

J 

and inxeger. 

This can alternatively be stated as, 

Itin. C X 

s.t. 

A X > b 

X >0 and integer, 

where C, b are a.ntegral vectors and A is an m x n matrix with 
all entries either zero or one. Without loss of generality 
w£ assme that b > 0, C >_ 0, and that A has atleast one 1 in 
each row. Each b^ may be imagined to represent the number of 
\‘;orkers required during time period i, and each c the cost of 
assigning a worker to shift Each column of A may be imagined 
to represent a possible work shift. 

In this dissertation, we consider a special case of the 

above problem that is c = 1, W j = 1,..., n. The problem then 

J 

IS to determine smallest number of v/orkers to satisfy the demand 


schedule. 
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Mathematically the problem is» 

Minimize IX 

s • t • 

A X b 

X >0 and integer. 

Subsequently this problem is referred to as cycla^c^ 
staff scheduling^ (CSS) problem. 

1.4 Lit eratiQ? e^ Jley 1 ew = 

The cyclic staff scheduling problem has been studied by 
several researchers in the literature. Dantzig L7j formulated 
the general CSS problem as an ILP as given belov. 


Min. 

n 

Z 

1=1 

°1 

uv 

1 


s.t. 

n 

Z 

1=1 

a 

11 

J\. 

1 

>, b^, 1 = 1, . . . , m 


•V 

1 


0. 

integer, o = 1, . . . , 


where X is the number of workers in shift j, and C is the cost 
J J 

of assigning a worker to shift a^ =1 if shift j is present 
on i“th time period and a. =0 otherwiso. The b is the require- 

J_ J a. 

ment during time period i. 

Henderson [9j has considered a special case of Dantzig 's 
general ILP formulation for scheduling of telephone operators 
as given below. 
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n 

Min. B 

X 


H 

1! 

0 


s * "b ♦ 

n 

z 

a X > b . 

1 — m 

J=1 

13 3 1 


X 

3 

0, integer, 



Luce [12] has suggested an alternative approach for 
scheduling of telephone operators. He essentially modified the 
ILP formulation of nhe CSS problem suggested by Henderson by 
adding shortage and surplus tolerance variables. Luce's ILP 
problem is more difficult nhan solving Henderson's problem as 
more integer variables are needed for formulation. However, a 
constraint fixing the total number of operators can be added in 
Luce's model. The underlying assumption of Luce's model is that 
"cost" of deficit equals "cost" of surplus. Penalizing under- 
staffing and over-staffing equally is indeed a little peculiar 
and unrealistic. 

The cyclic staff scheduling problem is more easily 
solvable when resources are continuously^ available. Baker [3] 
has suggested a method which needs only hand computations to 
solve the i5--7) cyclic staffing problem optimally. In this 
problem, full-time workers are available for five consecutive 
days in a one week cycle with two consecutive days off. 

Baker [4] has suggested another simple algorithm for the 
problem of allocating staff in order to meet demand over a 
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repeating seven-day cycle when both full time and part-time 
workers are there. The formulation assumes that; full-time 
employees are entitled to two consecutive days off and that 
the objective is to meet demand for staff with a minimumu use 
of part-time employees. 

For solving the general cyclic staffing problem with 
continuously available resources, Bartholdi et al [6] suggested 
an optimal technique based on linear programming round-off. In 
this method, the given CSS problem is first solved, as an LP by 
relaxing the integrality restriction on all the variables. The 
continuous valued optimum solution obtained is then roundea-off 
in a special wa 3 )-, which guarantees the feasibility?- and optimality 
of the rounded-off solution to the original problem. 

In practice, resources may not be available continuously. 
This may be due to lunch or relief breaks for the workers. The 
CSS problem with intermittently available resources is more 
difficult to solve. Current applications resolve the dilemma 
of problem size in a variety of ways. One possibility is to 
solve the sCSS problem as an LP and then to use some heuristic 
procedure to round-off the resulting non- integer solution. 

The schediiLiiig of intermittently available resources has 
been studied by Glover and Mulvey [s], Segal [iS], Henderson 
and Berry [lO], Laportc [ll], Shepardson and Marsten l19], 

Bartholdi [ 5 ], Ondhensden and Wen-Jang [I 6 J, Willis and Wilson [2l] 
Habert and Watts [l3j, and Morris and Showalter [i4] . 
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Segal [is] considered the problem o£ scheduling telephone 
operators where operators v/ork in shifts that arc specified by start 
and end times. Each tour contains from one to three relief periods 
during which the operators are not working, and a particular sche- 
dule of relief periods within a tour is called a trick. He 
proposes a method for determining the assignment of operators to 
each trick, it utilizes the out-of-Kilter algorithm and tvro 
network- flov/ models. The problem is formulated as an IIP. 
n 

Min, I C X 
0=1 0 0 

s.t. 


n 

I 

0=1 


^ ■y 

'10 0 


X 


0 




>, V i = 1, . . . , m 

0, integer, ¥ o = !»•••» 


n. 


He considers this problem unsolvable with the available 
integer programming techniques at that time. He follows an 
approach based on network-flow for formulation in which integer 
solutions arc guaranteed. Initially breaks and relief periods 
are ignored and the simplified problem is solved using out of Kiltei 
algorithm. He then brings in relief and break periods and uses 
some heuristic procedure to meet shortage at break and relief 
periods. 


Henderson and Berry [lO] gave some heuristic methods for 
scheduling of telephone traffic exchange operators to meet demand 
that vary over a 24-hour operating period. Two types of heuristics 
are described? 
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(1) for determining the work shift types to be considered 
in preparing an operator shift schedule and, 

(2) for constructing an operator shift schedule from a given 
set of work shift types. 

They give tv;o heuristics for determining the set of shift 
types. Once the working subset (shifts) has been defined, the 
given CSS problem is solved as a relaxed LP and all the variables 
are rounded-up and then tv/o heuristics are applied to remove the 
excess operators. 

Bartholdi [5] has suggested a guaranteed accuracy round- 
off algorithm for solving the CSS problem with intermittently 
available resources. We describe his method someT,;hat in detail 
in a separate section, 

Oudheusden and Wen-Jenq C16] consider the following 
problem, given operator requirements during different time periods 
of the day, a certain nmber of operators are assigned to the 
particular work schedules in order to meet, as well as possible, 
the varying demand for operators. If the number of operators 
for a day are fixed at p, the following model is suggested. 

Max. Y 

s.t. 

n 

( E a. X - b.) > Y, 1 = 1,..., m 

la 0 1 - ’ » » 

n 

E X = p 
1=1 ^ 

X 0 , integer, 0 = 1 ,..., n, 

J 
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Y is an integer variable which can be positive, zero or negative. 

This IS a maxmin problem as it maximizes the smallest 
difference between supply and demand in any period. When 
(optimum objective function value), takes the value -1, at least 
one time interval v/ill be understaffed by one operator. When 
equals 0, all demands are met and overstaffing in all 
periods is not possible, '/flien equals + 1, all operator 

requirements are met with atleast one extra worker. Results of 
one realistic study, applied to Taiwan Telecommunication adminis- 
trative Bureau are presented. 

Arthur and Ravindran [ 2 ] consider a problem of scheduling 
nurses under multiple object3.ves. They make use of a goal 
programming model which allows for consideration of multiple 
conflicting objectives inherent in scheduling a nursing staff. 

Willis and V/ilson [2lj present a case study of scheduling 
telephone betting oporators. They demonstrate the application of 
linear programming and network flow techniques for determining 
optimal shift sizes for telephone betting operators. I'labert and 
Watts [l3l present a simulation analysis of tour shift construction 
procedures. The^r provide some managerial suggestions on the 
issues of effectively staffing personnel on a weekly basis, 

Morris and Showalter [l4] give an exact and a heuristic 
algorithms for solving the CSS problem with intermittently 
available resources. First they consider the problem of 
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scheduling workers on a day, where the workers vrork for 8 hours 
with a one hour break after 4 hoijrs of v/ork in 24“hotrr cycle. 
For solving this (24 x 24) problem optimally, they suggest a 
mixed cutting plane and branch and bound algorithm. They give 
some heuristic approaches for solving the weekly scheduling 
problem of (l68 x 168) size. 

1.4.1 FarthpldiJ s^ Method . 

Bartholdi [5] has suggested a heuristic algorithm to 
find a feasible solution for the CSS problem when there are 
intermittent breaks in the shift. In this section, we describe 
his method in detail. 

Lot A be an arbitrary m x n 0-1 matrix whose i-th rovr is 
and \/hoso j-th column is A, Suppose each columin of A has 
atloast K ones but no row has more than q blocks of consecutive 
ones (v/hon considered circularity). A vector whose every entry 
is (q “ l) is denoted by (q-l). 

The algorithmic steps involved in his mechod are given 

below. 

Stpp^l^; Solve the continuous-valued linear program LP(b+q-l)i 

I4ii. 1 k 

s. t. 

Ax b + q — 1 

X > 0 

Let the solution bo Xj^p (b+q-l) = (x^, Xp,,.,, 
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Construct the (integer- valued) 

— ( r ^ ^ ^ ^1 **"^2 ^ 


heuristic solution, 

n n-1 

[2x1-; 2x1). 

h=l ^ h=l ^ 


Bartholdi has proved that the solution obtained after 
rounding-off as described above is feasible to the following 
problem. 


Min 1 X 

s.t. 

A X > b 

X 0, integer 


He also obtains an upper bound on the possible deviation 
of his heuristic solution (Z^^) from the integer optimum solution 
Zjp to the staffing problem. This boLind is, 

\ “ 2jp < i I Cq - 1) i. 

It might be noted that this bound is not satisfactory 
for higher values of q. 

1,5 Outl the^ jhp s is^ s 

A brief chapter by chapter outline of the thesis is given 
in this section. 


In Chapter II, we develop a heuristic algorithm for the 
CSS problem. In xhis method, the given staffing problem is 
solved as an ordinary LP by relaxing integrality restriction 
on all variables. An irnoroved version of the special rounding- 
off procedure suggested by Bartholdi et al l6] is applied. If no 
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constraint is violated with the rounded~off solution, this 
solution IS optimum. Otherv/ise, v;e find out the unsatisfied 
demand in each time period. This demand is taken as the require- 
ment of the new CSS problem. This is repeated until there is no 
unsatisfied demand for every time period. The solution obtained 
Using this algorithm is found to be very close to the optimum 
solution. A numerical example is presented to illustrate the 
algorithm. 

In Chapter III, vre develop an adaptation of the general 
integer programming branch and bound algorithm and a dual 
fractional cutting plane algorithm for the CSS problem. Some 
of the important features of this algorithm are. 

(i) Only one A-matrix of fixed size is needed to be stored 
throughout the problem, 

\,ia ) Reoptimization is greatl3^ used to reduce the computations 
substantially, 

(ill) To unloose lower (upper) bound t^qie of constraints on an3?- 
variable, sensitivity analysis of the lovrer (upper) bound is 
performed, 

(iv) The heuristic developed in Chapter II is applied at every 
node of the branch and bound tree to obtain a good feasible 
solution, 

(v) The variable selected for branching is such that the 
pruning is more effective. 
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Detailed discussion of all these aspects is presented, A 
numerical example is presented to clarify the above steps. 

¥e also develop a cutting plane algorithm for finding an 
optimum solution to the CSS problem. Dual fractional cutting 
plane algorithm is used. Two types of cuts are introduced. It 
IS known that ixi an integer optimum solution to the CSS problem, 
sum of all variables must be an integer, so whenever the objective 
value IS found to be fractional, we apply the cut that the sum of 
all original variables must be greater than or equal to the least 
integer greater than the objective function value. If the objec- 
tive function value is fractional one, wo then apply Gomory 's 
fractional cut. Detailed discussion about the strateg3^ used 
for finding source variable for cut generation is also given. 

In Chapter IV, we present the computational performance 
of Branch and Bound algorithm, LP-based heuristic, cutting plane 
algorithm and Bartholdi's method. The results are found to be 
quite encouraging. The solutions obtained using LP-based heuris- 
tic were found to be very close to the optimum solution. The 
number of LPb solved in any problem were never more than 2. The 
Branch and Bound algorithm is able to solve problems of size 
(96 X 96) optimally in a fev; minutes of CPU time. The perfor- 
mance of cutting plane algorithm was observed to be erratic. It 
could not solve problems of size more than 1.48 x 48), Computer 
programs for all these algorithms were written in FORTRAN- IV and 
implemented on DEC-IO9O Computer System. 



CHAPTER II 


LP - BASED HEURISTIC AIDORITHI'I 
2 . 3- jjrbrodi^ti pri i 

In this chapter, a new heuristic algorithm for solving 
the CSS problem is developed. The algorithm is called IP-b a aed 
heuristic algorithm as a sequence of IP's are solved to get a 
good feasible solution to the CSS problem. ¥e first describe 
the intuitive ideas that led to its development, A step\/ise 
description of the algorithm is then presented, A numerical 
example is solved to illustrate the algorithm. 

2,2 Development of the Heuristic Algorithm I 

Bartholdi [5] developed a heuristic algorithm for solving 
the CSS problem. He gave a priori bound for the maximum possible 
deviation of his heuristic solution (Zj^) from the integer optimum 
solution (^jp) I'or uhe CSS problem. The bound is 

“ '^IP ^ ^ k ~ l)"' • 

His algorithm does not give a good solution for the CSS problem 
when there are large nunber of intermittent breaks (i.e,, 
Bartholdi's q is large) and/or when n/k ratio is large. For 
example, if \/e consider a (96 x 96 ) problem with k = 30 and 
q = 6, Bartholdi's bound is i (6~l)l = I 6 units. Further, 

it was observed by us computationally that his algorithm 
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invariably requires as manji^ additional imits as given b 3 '' his 
bound because he is adding (q-l) additional units to the actual 
demand in each time period in the initial stage to obtain a 
feasible solution to the original CSS problem. «/e observed 
that there is no need to add (q-l) units to the requirement in 
each time period in the initial stage itself v;ithout examining 
whether it is really needed or not. Our strategy is to provide 
the additional number of v/orkers on those days where the need 
arises after appl^rmg the special rounding-off to the first LP 
solution obtained. With this intuitive idea in mind, we develop 
the LP“based heuristic algorithm, for solving the CSS problem. 

Familiarity \uth the following special v/ay of rounding- 
off IS needed before going through the formal description of the 
algorithm. Given n continuous valued variables say X 25 ...,x ^9 
the^/" can be rounded-off in the following special wa}/' suggested by 
Bartholdi et al [6]. 

_ n n-1 

Y - (ix. 1, [x-j-i-Xpl - [x. 1 , , . . , r s X ] - r Y X 1) 

J=1 1=1 ^ 

(?.2.1j 

Bartholdi has shorn that il a CSS problem is solved as an 
ordinary LP and the continuous n values obtained are rounded-off 
in the above fashion, there ■'nill not bo more than (q-l) units 
of unsatisfied demand on any of the time oeriods. We generalize 
this result to the n integer solutions obtained by the following 
round-offs. 
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To find-out whether a given solution is feasible to the 

given CSS problem or not we find out the unsatisfied demand for 

each time period with thct solution. If the unsatisfied demand 

on all the time periods is zero, the jiiven solution is feasible 

to the CSS problem. Let -- '^2^* •••> Y-^l denote the 

rounded“Off solution vector corresponding to the h-th rounded-off 

h 

solution, where denotes the rounded -off value of the variable 
in the h“th solution vector. Let 1''^ = iu^^, 






- 


v/here, = ^ 0 -n’ 0 > 
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h. 

denote the vector of unsatisfied demands -v/iien Y is taken as the 

*h 

solution for the CSS problem, denotes the unsatisfied demand 

■U 

on i“th time period when Y'^'^ is taken as the solution for the CSS 

h 

problem. Each u^ is obtained as given belowi 

T1 

u ^ - maximum { 0, b - Z a ,y (2,2,2) 

Vi=: 1,,.., m, ¥h = 1,...5 n 

Now, we explain how the LP-based heuristic proceeds to 

get a good feasible solution to the CSS problem. V/e first ignore 

the integrality/- restriction on all the variables and solve an 

LP to get continuous valued optimu^i solution to the CSS problem. 

Now, n- rounded-off solutions Y^, Y^, Y^ and the corresponding 

12 n 

unsatisfied demand vectors d , U , ..., U are obtained. Total 

h h 

unsatisfied demand for solution Y^ is given by Z u , We 

1=1 ^ 

calculate the total unsatisfied demand for each of the n solutions 
obtained and select the one I'hich gives the minimum amount of 
total unsatisfied demand. Let this solution be Y^. If the 
total unsatisfied demand for this solution is zero, we have a 
feasible solution to the CSS problem and stop. Othenvise, a new 
LP is formulated by taking the demand vector for the given CSS 
problem as uP, Wo now solve this LP and obtain n ro-unded-off 
solutions and select the one with miniium amount of total unsatis- 
fied demand. The values of all y P' s, of the best solution 

J 

obtained here are added to the corresponding rounded-off y P 's 

tj 
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obtained upto the previous LP, This process is repeated until 
the total unsatisfied demand becomes zero. 

A formal statement of the algorithm is given in the next 
section. 


2.3 siemprt_ of ^ tb6_ 

Step^ m Set X* = 0 . 

Step^_l^. Let X be the solution of the following continuous valued 
variable CSS problem. 


Min. 1 X 

s • t . 

A X > b 
X > 0 

*1 ^ n. 

Step^ Obtain the n round-off inreger solutions Y , Y , . . . , Y 
using (2. 2.1.1) to (2.2.1,n). Compute tho vectors of 
unsatisfied demands using (2.2.2) correspoiiding 

to Y^j Y^, Y^ respectively. Let be the solution 


for i/iiich 

m m 

Z u^ < 2 n. ¥3=1,..., n. 

i=l ^ " 1=1 ^ 

Set X * = -f yP . 

m 

If Z u.^ = 0 go to Step 3 j otherwise set b = and 
1=1 ^ 

go to Step 1. 

S_tep_ _ 3 s. The vector X'^ is the heuristic optimum solution of the 


given CSS problem. STOP. 
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It has been computationally observed that deleting all 
those constraints for which b^ is zero while solving the LP in 
Step 1 results in significant reduction in computational time. 
2,4 Nmer 1C al^ * 

Me now solve a numerical example to iHlustrate various 

steps involved in the LP -based heuristic algorithm. Consider the 

follov/ing hyoothetical CSS problem. 

12 

Min. 2 X 
1=1 

s.t . 




'l 

0 

0 

0 

0 

1 

0 

1 

1 

0 

1 

I'l 1 
\ 

^1 


123 

1 

1 

0 

0 

0 

0 

1 

0 

1 

1 

0 

1 


X2 


123 

1 

1 

1 

0 

0 

0 

0 

1 

0 

1 

1 

0 


X- 

p 


122 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

1 

1 


^4 


125 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

1 


^5 


123 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 


^6 


124 

0 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 

X 


x? 1 


124 

1 

0 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 


^8 


122 

0 

1 

0 

1 

1 

JL 

0 

1 

1 

1 

0 

0 

0 




123 

0 

0 

1 

0 

1 

1 

0 

1 

1 

1 

0 

0 


^10 


124 

0 

0 

0 

1 

0 

1 

1 

0 

1 

1 

1 

0 i 


^11 


124 

*0 

1 

0 

0 

0 

1 

0 

1 

1 

0 

1 

1 

1 


^"12 


123 

\ 














^ 1 

1 


X >, 0 , integer, V 3 - 1 ,..., 12 . 

cJ 
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This problem is first solved by relaxing the integrality 

restriction on all x 's as an ordinary LP. The optimal 

J 

continuous valued variable vector is» 

X = (39.80, 0.80, 40.60, 0.80, 41.40, 0.40, 40.40, 

0.00, 41.60, 0.00, 40.80, 0.40) 

Step^ ps The n rounded-off solutions ares 


y1 


(40, 

1, 

Al, 

0, 

42, 

0, 

41, 

0, 

4ls 

0, 

41, 

0) 

y2 

= 

(39, 

1. 

41, 

1, 

41, 

0, 

4l, 

0, 

41, 

0, 

41, 

1) 

y5 


(40, 

0. 

41, 

1, 

4l, 

1, 

40, 

0, 

42, 

0, 

40, 

1) 

y^ 


(40, 

1, 

40, 

1, 

42, 

0, 

AO, 

0, 

42, 

0, 

41, 

0) 

y5 

r= 

(40, 


4l, 

0, 

42, 

0, 

4l, 

0, 

41, 

0, 

41, 

0) 

y6 

= 

(AO, 

1, 

40, 

1, 

41, 

1, 

40, 

0, 

42, 

0, 

41, 

0) 

y7 


(39, 

1, 

41, 

1, 

41, 

0, 

41, 

0, 

41, 

0, 

41, 

1) 

yS 

= 

(40, 

1, 

40, 

1, 

42, 

0, 

40, 

0, 

42, 

0, 

41, 

0) 

y9 

=: 

(40, 

1, 

40, 

1, 

42, 

0, 

40, 

0, 

42, 

0, 

41, 

0) 

ylO 


(39, 

1, 

41, 

1, 

41, 

0, 

41, 

0, 

41, 

0, 

41, 

1) 

yll 


(39, 

1, 


1, 

41, 

^9 

41, 

0, 

41, 

0, 

41, 

1) 

y12 

rr 

(AO, 

0, 

Al, 

1, 

41, 

1, 

40, 

0, 

42, 

0, 

40, 

1) 


The unsatisfied demand vectors and the total 
unsatisfied demand for each of the above 12 solutions 
are given in Table 3.1. 
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day 


U. 


1 


U. 


u. 


u 


4 


U, 


U, 


U 


7 


u. 


8 


U. 


“lO 


u 


12 


123456789 10 11 12 T.U.S. 

. _^. ... . 


1001000000 1 0 

1000010001 0 0 

0010010000 O 1 

0000001100 0 0 


3 

3 

3 

2 


1001000000 1 0 

OOOOIOIOOO 0 1 

1000010001 0 0 

0000001100 0 0 


3 

3 

3 

2 


0000001100 0 0 

1 000010001 0 0 

1000010001 0 0 

OOlOOlOOOO 0 1 


2 

3 

3 

3 


Table 3.1 

We find tho.t the total unsatisfied demand (T.U.S, demand) 
corresponding to thesolution is minimum. Therefore, 

= (40, 1, Z^O, 1, 42, 0, 40, 0, 42, 0, 4l, 0) 

Step_35‘ The new demand vector for the CSS problem is 

b = ( 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0). 

Solve the following LP, with this demand vector. 
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12 

Mxn, I 2 : 
0=1 ^ 

s.t. 


.-2 -3 ' -5 • -5 ^ ^ 

‘h X^, Xg -!- X^ -- Xg 1 

^-1 > 0 * 

O 

After solving this LP, we get, 

X = (0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0) 


Since this is an integer solution, if we round it off 
in 12 V/83/-S, we get the same solution. The total unsa- 
tisfied demand is zero v/ith this solution. Thus = X. 
Therefore, 

X* = -f- yP 

= (40, 1, 41, 1, 42, 0, ^lO, 0, 42, 0, 41, O) 


This is the final solution to the given CSS problem. 
In this esse, this is also an optimum solution. 



CHAPTER III 


EXACT AIC-0RITM3 


3 . 1 jffltrodjj.pt ipn . 

In this chapter, we develop tv/o exact algorithms for 
solving the CSS problem. In Section 3.3, we develop a branch 
and bound algorithm for solving the CSS problem optimally. This 
algorithm is an adaptation of the general integer programming 
branch and bound algorithm. Some of the important features of 
the proposed algorithm ares (l) Sensitivit^A analysis is made to 
impose lower or upper bound type of constraint on any variable 
to signif 1 C antler reduce the storage and computational time 
requirements, (2) The LP-based heuristic algorithm developed 
in Chapter II is applied at every candidate problem to get a good 
feasible solution. ¥e first explain in detail, different strate- 
gies that are follov/ed in the algorithm. The flovz-chart for theps: 
proposed branch and bound algorithm is presented. A step-by-step 
procedure, to perform sensitivitj^ analysis for imfflosing lower 
(upper) bound constraint on some variable, is then given. A nume- 
rical example is also solved. In Section 3.4, v/e develop a dual 
fractional cutting plane algorithm for solving the CSS problem. 

We introduce two types of cuts in this algorithm. One cut is 
based on the objective function value and the other one is 
Gomory's fractional cut. 
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3 . 2 Nota tions » 

To Simplify discussion in this chapter, we make use 
of the follov^’ing addicional symbols and notations. 

Z optimum objective function value of the continuous 

valued CSS problem. 

Z^ lower bound of the solution for the CSS problem. 

UB upper bound of the solution for the CSS problem. 

Z" incumbent objective function value. 

X incumbent solution. 

continuous valued solution of candidate problem P^. 

Z^ optimum objective function value of the continuous 

valued candidate problem P^. 

Zj^ objective function value of the integer feasible 

solution obtained when the LP-based heuristic algorithm 

is applied by taking as the first LP solution, 
k 

Xj Integer feasible solution obtained when the LP-based 

heuristic is applied at candidate problem P^. 

L array containing lovrer bounds on variables 

x^, 1 = 1, . . . , m+n. 

U array containing upper bounds on variables 

x^, 1 = 1, , . . , m-rn. 

B {bl, b2, ..., bm} set of all basic variables, v/here 

bi denotes the index of the basic variable corres- 
ponding to the i-th row. 
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3 . 3 Bra nch._a ji_d._ Bo ' 

The basic idea in this method is to investigate all the 
feasible solutions by partitioning the solution set into subsets 
by a process, knowi as branching. A lower bound is calculated 
for the solutions with in each subset. If the lower bound is 
smaller than the objective function value of the incumbent 
solution and the subset with the lower I'ound is feasible, further 
branching is stopped from this subset. This is known as fathom^ 
ing. After each partitioning, those subsets with a bound that 
exceeds that of a kno^yvn feasible solution are excluded from 
further partitionings. This is knoi'n as prming. The partitioning 
continues until a feasible solution is iound such that its value 
IS not greater than the bound for any subset. The algorithm 
ternn nates when all the subsets are either pruned or fathomed. 

We now describe various strategies of the branch and bound algo- 
rithm for the CSS problem considered in this thesis. 

A branch and bound algorithm is cnaracterized by its 
(i) bounding strategy, (li) fatiioming strateg 3 '-, (iii) pruning 
strategy, (iv) branching strategy, and (v) searcii strategy. 

These strategies for the suggested branch and bound algorithm 
for the CSS problem are discussed in the follo^'^ing sub-sections. 

3.3.1 Bp undi ng^ S t r at egy ' 

First, we solve the given CSS problem by relaxing the 
integrality'' restriction on all variables. The lower bound for 
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the given CSS problem is the smallest integer value greater than 
or equal to the continuous valued obooctive function. 

LB = Z = [Z ]. 

The lov;-er bound for candidate problem = I Z^l . 
Initially we do not have any incumbent, therefore, 

Z = oo . 

5.3.2 Fathoming^ Otj'at e_^j 

Wc fathom a candidate problem if Z^^ satisfies the 
following condition. 

r Z^l = Zj^^ 

3.3.3 Prymang^ S;brat egy 4 

We prune a candidate problem P^, if either of the follow- 
ing conditions are satisfied 

[Z^l Z , or 

candidate problem P^ is infeasible. 

3.3.^ PSY * 

1 /' 

If a candidate problem P^ is neither fathomed nor pruned, 
it is used for subsequent branching. To perform branching at a 
candidate problem P^, we can follov; different strategies. The 
variable selected for branching shoudd bo such that it results in 
the enumeration of minimum number of nodes for finding the opti- 
mum solution. Many strategies \/ere computationally tried by us. 
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The follov;ing strategj^ produced consistently better results 
than the other strategies. 

As observed earlier in Section 2.2, by rounding-off the 

solution obtained at candidate problem using (2. 2,1,1) to 

(2,2,l.n) we get n integer solutions. Let be the solution 

with minimum amount of total unsatisfied demand and IJP be the 

unsatisfied demand vector corresponding to this solution. If 
m 

I u P is zero we fathom this candidate problem otherwise we 
1=1 ^ 

select a variable x for further branching. The strategy used 

Id 

for selecting x_ is explained below. 

o 

have available vath us. Now, find out, F , the 

J 

number of days on which shift i is active during the days for 
which demand is not satisfied. For each of the basic variables 

we compute, f = ix 1 - x , ¥c select a variable x such that 

J J J s 

F X f > F X f , V J = 1,..,, n. It is observed that 

o S J J 

selecting a variable suen as x for branching brings in maximum 
amount of change in the lo\fer bound in most of the cases. Now, 
we give a step by step procedure to identify/- branching variable 

^s- 

Step_ 05 Seb E = F = f = 0, V j = 1,..., n 
* ** J J J 

Step, 1, For all i = 1,..., m if (u^-^ >0) then 

for all 0 = 1 to n such that a =1, set F = F +1, 

for all 3 = 1 to n such that j s B, set f = Fx 1 — x . 

J J J 

For all 0 = 1 to n such that 3 £ B, set E = F x f , 

U J J 
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Find out a variable x such that 

E = maximun (E^f 
^ l<d<n ^ 

Select the variable x as the source variable for further 

o 

branching * 

It IS computationally observed that imposing upper bound 
type of constraint first results in the enumeration of smaller 
number of nodes for most of the problems. Hence, this strategy 
has been implemented. 

3.3.5 S t r at ? 

As far as searching strategy is concerned, depth first 
search is employed because of the ease in implem.entation and 
lesser storage requirements. 

We give the flov; chart of the proposed branch and bound 
algorithm in Fig. 3.1. 

3.3*6 Sensitiyitj^ s 

In this section, we present a step by step procedure to 

perform sensitivity analysis by implicitly considering lower 

bound or upper bound type of constraint on some variable 3t . Let 

B denote the set of all basic variables ab any stage of the 

problem. Corresponding to this basis, the relative cost Hcoeffi- 

n 

cients of the variables x^, j =1,..., (m-Ki) are c_ = c - Z 

j ^ ^ X=1 

Set C = {gs c 0} and H = {j. c < 0} where G and H respectively 
denote the variables at their lov/er bound and at their upper bounds 



S5 s» 00 I k « 1 , T 




k « k-^ t 

1 


Solvt candidate probl^ P 
Min. 1 X 

A X 3j b 
X 2 0 

I.(r)< U(r), r- 1,..,I 


I - TUI 



p fMuiibi. 

1 y.« 


k k* If 

I»«t X b« tha optimum solution to P with 2 ms 
ths objsetlvs tunetlon Talus. If k» 1 than 

r . PI. 

AppXX ths LP^bastd hsurlstls slgoritlin to 

obtain an Intsi^sr fsaslbls solution 
k k 

Xj with Zj as ths objsotiTS funetlon valus. 



V< 2 

^ k 


2 H Zj 
# k 

X « Xj 


z"* > 7 


flnd-out th« bruiobtnu varlkbl. 
t - T+1, 8(T) - •) i(r) ■R ) U{T) . 0 

T <• T+1} 3(T) - ■) L't) - 0; U(T) . \j'» J| 
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Let (B, G, H) be an optimal basis structure of the CSS 
problem. Let A = la } = A. Further, let be the 

amount of change required in the current value of the variables 
Xg to satisfy the bound imposed on it. We first try to make 
the variable x non-basic, if it is a basic variable. If x_ is 
a basic variable and no variable is found to be el.igible to 
enter the basis, the problem is infeasible. If s e B and there 
is some variable t eligible to enter the basis, we make the 
variable x non-basic by performing one dual simplex iteration 
C 20 J, Now, we find out the maximum amount of required change 
that can be brought in the value of the variable x without 
changing the present basis structure. We denote this possible 
change by/ . If is greater than the required change S > ’W’e 
bring~in5 amount of change in the value of the variable x and 
update the b column for the corresponding change and stop. If 
^ is less than 5^, orly A amount of change can be made in the 
value of the variable x ’■’Uthout changing the present basis. We 
bring A amount of charge in the value of x and update b column 
for the corresponding change brought in x , In this case, the 

0 

bound restriction on the variable x^ is not yet satisfied because 
one of the basic variables reaches eitner its lower bound or its 
upper bound before the required change in the value of x can be 
brought in. The dual simplex iteration is again performed to 
replace the basic variable, which has reached its lower or 
upper bound recently, by an appropriate non-basic variable. The 



55 


basis structure is updated and tins process is repeated until 
either = 0 or it is found that the dual simplex iteration can 
not be performed indicating infeasibility of the problem. A 
formal description of the procedure to perform sensitivity 
analysis is given belov/, 

Vfe consider the case when some lox/er bound restriction 

IS imposed on variable x , V/o give the necessary changes 

s 

reqnred to be done if upper bound type of restriction is to be 
imposed on x„. 

o 

step, p; Let B denote the set of all basic variables, 

G = Ijs c >, O} , and H = {jc c" s 0}^ where, 

J J 

m 

c' = c - 2 a , ¥ 3 = 1 ,..., n. Let L and U be two 

^ 1=1 

arrays respectively containing the tightest lov/er bound 
and upper bound restrictions on each variable 3, 

¥ 3 = 1, . . m+n. 

Set L(s) = X . 

If s / B go to Step 2. 

Let r be the rov/ corresponding to x in A, 

Step ^13 Tho variable leaves the basis at its lower bound 

or upper bound. If x-^^ is equal to L(br) go to l(a), 
else go to l(b). 

l(a)i Let K = {3S 3 e G and a <0, or 3 £ K and 

J 

a >0}. If K ms empty, go to Stop 3, otherwise 
r 3 

identify a non~basic variable x^ for which 
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-c./a , = min. {-c /a' }, and go to l(c). 

^ oeK ^ 

1(13)5 Let K = (35 3 e G and a > 0 , or 3 e H and 

■1 J 

a < 0} . If K IS empty, go to Step 3 , otherwise 

I J 

identify a non-basic variable x.j. for which 

c./a ^ = min. (c /a }, and go to l(c). 

3eK ^ 

1(0)5 The variable leaves the basis and x^ enters 
the basis. Update the matrix (except the column 
corresponding to by performing one dual simplex 
iteration. Go to Step 2. 

St ep 2 5 Update all x e B as follows S 

J 

Compute S' = rL(sTl - L(s) 

If B as zero go to Step 4 , otherwise 
compute. 



' (S-^ - L(bi))/a.g , If > 0 

] “ > hs = ° 

\^-(U(bi) - t )/aL , if a < 0 


Let - min. { A , S} 
l£i<m 

Now update = x^^ - A 

L(s) = L(s) + A 


If 5 = ^ go to step 4 otherwise identify the row r 
corresponding to the minimum and go to Step 1. 
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St ep The problem is infeasible. STOP. 

'I'he solution X is optiraum and STOP. 


For performing sensitivity analysis to impose upper 
bound type of constraint on variable x , we make the following 
changes in Step 0 and Step 2. 


Step^ 


Set U(s) = X instead of L(s) = x . 

s s 

Update all x e B as follows. 

J 

Compute ^ = U(s) “ lU( s ) J 


If ^ is zero go to Step 4, otnerwise, 
compute, 


Ai = 




{u(bO - b'p/a;^ 


V 


■(b. - L(bi))/a 


IS' 


if a 


IS 


if a = 

IS 

if a < 0 
X s 


> 0 
= 0 


Set A = min } 


KKm 


Now update x^. 


X-, , + a 
bi IS 


A 

u(s) = U(s) - A 


s s 


If go "to step 4^ otherwise identify the rov/ r 

corresponding to the minimum and go to otep 1. 

If a node in the branching tree is either fathomed or 
pruned v/e have to get back to its parent, for this standard 
bounded linear programming techniques [lb] can be used. 
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The algorithm is coded in FORTfliN-XV on DEC-1090 computer 
system and the listings are enclosed at the end. Extensive 
computational investigation is done. The performance of the 
algorithm is quite encouraging and is discussed in Chapter TV. 

We nov/ take up a numerical example to illustrate the method. 

3.3.7 Numerical^ ° 

We take up the same example that is considered in 
Chapter II. The branching tree enumerated is given in Fig. 3.2, 
Nodes are numbered by the order in which they are enumerated. 

The branching variable and the bound imposed on it is shown at 
each nodo. The lower bound obtained at each node is written in 
a rectangle at that node. Initially, the given CSS problem is 
solved as an ordinary LP, the continuous valued objective func- 
tion value obtained is = 247.00 [z^l = [247.00] = 247 is the 
lower bound for the CSS problem. When LP-based heuristic is 
applied at node 1, we get a feasible solution with Zj = 248. 
Therefore, Z = 248. The algorithm terminated after enumerating 
25 nodes. In this problem the initial solution obtained at 
node 1 itself is the optimum solution. 

X = (40, 1, 41 , 1 , 42, 0, 40, 0, 42, 0, 4l, O) . 
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3*4 Cuttiiag Plane Algorithm s 
3.4,1 Introduction i 


In this section, we develop a dual fractional cutting 
plane algorithm for solving the CSS problem. The general 
intent of cutting plane algorithms is to deduce supplementary 
inequalities from the integrality and constraint requirements 
and to append to the constraint set so as to eventually produce 
a linear program whose optimal solution is integer In the 
integer constrained variables. Gomory [17] developed the 
first cutting plane algorithm applicable to any Integer Program. 
Gomory [l7] produced a second cutting plane algorithm for the 
integer program which requires only additions and subtractions 
in computation (an " all- integer " technique) . All these methods 
maintain linear programs vdiich are dual feasible, and are there- 
fore often classified as dual cutting plane algorithms. 


Glover and Young developed cutting plane algorithms for 
the Integer Program which maintain linear programs that are 
primal feasible. Since primal feasible integer solutions are 
successively produced, the technique is referred to as a 
primal cutting plane algorithm. 


In this dissertation, we consider a cutting plane algo- 
rithm for the CSS problem which utilizes the lexicographic dual 
simplex method and allov/s fractional numbers in computation - 
hence the dual fractional reference. We take up the dual 
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fractional cutting plane algorithm because of its faster 
convergence rate. 

5*4,2 Pt. Algorithm, ° 

To avoid cycling and to ensure convergence we make use 
of lexicographic dual simplex method in the cubting plane 
algorithm. 

Consider the problem. 

Max. C X = Z 

s . t ♦ 

^ ^ ^ integer 

This problem is represented in the following way in an 
inverse tableau to make use of lexicographic dual sinpilex method, 



1 

-*x^ 

, . . . 

^ « 

“^n 



^0 = 

^oo 

%1 

• # * 


9 0 » • 

%n : 



Xi = 

0 

-1 


0 

» * « • 

0 I 

> 

0 

• 

• 







# 

• 


X . = 

0 

0 

0 

• » « 

“1 

• • • « 

0 

> 

0 

» 

• 

''n = 

0 

0 


0 

* « « * 

-1 

• 

• 

> 

1 

0 

^n+l = 

^n+1, 0 

n-’-l , 1 


^n+1,0 

« « • • 

®n+l,n 

> 

0 

• 

^n+m ~ 

cL 

n+m, 0 

‘^n+m,! 


3 

n+m, 2 

• « • « 

a , 

n+m,n 

4,^ 

• 

• 

0 
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where 




- Of 
00 * 

= “C , 

01 

¥ 

J - 

n 

, • = 

b 5 

¥ 

~ ^ 

m 

n+i, 0 





a ^ = 

a . 

¥ 

IL ~ 3- ^ ^ 

m, 

nH-1,3 




7 


* X ^ y 3^ 


In addition, we use the following notation: 


<J contains the n indices of the current non-hasic variables 
as they appear from left to right on the top most line of 
the tableau. Initially, J = {1,2,..., n} , 

'I(l) (o = l,,..,n) is the o-th element in J. 


Now, we explain how our cutting plane algorithm proceeds, 

We first solve the given CSS problem as an LP using lexicographic 
dual simplex method [l5] and then round--off the continuoiis 
valued solution obtained in n wa^rs using (2, 2,1.1) to (2.2,l.n) 
to obtain Y^. Let Y^ be the solution satisfying the 

follov7ing inequality: 

E u-P < E u.'^, ^3=1,..,, n. 

i=l ^ “ 1=1 ^ 


If UP =■ 0, then yP IS an optimum solution to the CSS 
problem and we stop. Otherv^ise, we impose a cut vdiich elimi- 
nates the current non-integer solution and is satisfied by all 
the integer feasible solutions to the given CSS problem. ¥e 
propose to impose two t 3 rpes of cuts. One cut is based on the 
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objective ftinction value (Z) of the optimum solution and the 

f 

other cut is Gomory s fractional cut. 


Since the CSS problem is a pure integer program, the 
objective function value of the optimum solution to this problem 
must be an integer. Therefore, \vhenever we observe [zl > Z, we 
append the following cut to inverse tableau. 


n 

k=l ^ 


‘Zl 


(3,4, 2,1) 


However, this cut can not be added directly in the above 
form. This problem can be tackled in the following way. Ini- 
tially, before solving the LP, we append one more additional 
n 

constraint Z x . > 0 in the (in+n^-l)-th row and then solve the 
0=1 ^ 

LP. Let slack variable corresponding to this 

constraint. We show that (3.4. 2, l) is equivalent to, 

n 


X 


rzl + Z a 


) 


( 3 . 4 . 2 . 2 ) 


m+n-i-1 “ “ ^Zi m+n+ 1,0 

J — J. 

The row associated v;ith x^ in the current inverse tableau 


leads to constraint, 
x^ = b 


V 


n 

+ Z a (“Xt / . \ ) 
vj ^ J(d) 


(3. 4.2. 3) 


0=1 


Therefore, from (3.4. 2. l) and (3. 4.2. 3) 


n 

Z b, 
k=l * 


n 

+ E 




( z a,) (-Xj/ > [z] 

0=1 0=1 ■ 


or, 
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n 

Z 

k=l 



n 

+ s 

3=1 


n 

( 2 
k=l 





therefore, 

n 

^ .f-i ^in+n+l,j ^**^1(3)^ ” ^m+n +1 

vJ “-J- 


[zl 


= ^ - I'Z'' * 

Therefore, whenever we have to impose a constraint like 
(3«4,2.1), v/e simply make q = (2 - fzl ), Now 

(m+n+l)-th constraint violates primal feasibility because 
(2“ - [Z] ) is always negative. Apply dual simplex method to 
make the tableau primal feasible. 

The second type of cut we impose is the standard Gomory'^s 
fractional cut, ¥e impose this cut when 2 is an integer, but 
none of the n rounded-off solutions are feasible to the given 
CSS problem. If x^ is the source variable the following 
constraint is appended at the bottom of the inverse tableau. 

=^(nn-n+l)rt " -^vo * J, <3.4.2.4) 

where, is Gomory^s slack variable associated with 

k-th fractional cut added. 


F - a - [a ^ 3 = 0,1,,.., n 

vj vj *- vj j J • > 

3.4.3 §oi>rp§, ^ 

Proper selection of source variable for #at generation 
plays very important role in getting a fast optimal solution. 
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Let be the source row. The hyperplane defining the cut is, 

ho ^J(o) = h- 

We may say the larger the value of the ratio 
the stronger the cut is Li?]. Hence, the objective when 
selecting a source rov; is to produce an inequality or h3rper-- 
plane in which these ratios are as large as possible. Since 
source rows will usually yield inequalities in which the 
ratios are large for some variables and small for others, it is 
not clear as to which row to select. However, we intuitively 
feel that the following strategy produces a deep cut. 




Stej3_ 2 . 


Set f^ = b4 - ^ i = 1,..., m-.-n 

and = 0. 

Arrange the f^^ s in decreasing order of their magni- 
tude and select the first t s. 

For each of the t f^'' s selected in Step 1 
n 


Compute 2 f and N. , the total number of non-zero 
a-i 


f ^ s in the i-th row. Set F = N 

ij ijLg 


Let F = max. {F.} 

^ l<i<m:-n ^ 

Select the variable as the source variable for cut 
generation. 


Wo have computationally observed that if the value of 
t is between 10 and 20, it woiold give better results. A formal 
statement of the algorithm is presented in the next section. 
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3.4*4 Stat ement^pl J:he_ i^^orithm : 

Solve the following LP using Lexico-dual simplex method 
after transforming the given CSS problem into the 
inverse tableau form. 

Min. 1 X 

s.t . 

A X > b 
n 

2 X > 0 

0=1 ^ 

X 0 

Let X be the optimum solution obtained. 

Use n round-offs (2,2,l.l) to (2.2,l.n) to get 
Y^, Y^. Select Y^ satisfying the following 

inequality^ 

Z u.P < Z u^, V j=l,...,n 

1=1 ^ 1=1 ^ 

If IjP = 0, go to Step 5, otherwise go to Step 3. 

Step^ If [z] = Z go to Step 4. Otherwise, introduce the 

cut (3, 4. 2. 2). Perfom necessary dual simplex itera- 
tions to get primal feasibility and go to Step 2. 

Find out the source row variable x^ using the routine 
described above and append the cut (3.4.2, 3) to the 
bottom of the inverse tableau. Perform the dual simplex 
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iterations to make the problem primal feasible. 

Go to Step 2. 

yP is the optimum solution to the CSS problem. STOP, 

The computational performance of the algorithm is 
discussed in Chapter 17. 



CHAPTER IV 


COr^PUTATIONAL INVESTIGATIONS 


4 . 1 ipduct s 

In this chapter, we present the •omputational performance 
of LP-based heuristic algorithm developed in Chapter II, the 
branch and bound algorithm and the dual fractional cutting plane 
algorithm developed in Chapter III and Bartholdi '^s guaranteed 
accuracy round-off algorithm. 

The computational experiments require a large number of 
problems. The details of the problems considered for computa- 
tional study arc given in Section 4,2. In Section 4.3, we 
discuss the various considerations such as accuracy, CPU time 
for analysing the three algorithms developed by us and Bartholdi 
algorithm. Finally, in Section 4,4, the concluding remarks are 
given, 

4.2 ^j3bl C ons idered, f or Expe rimentations 

Wo have considered five different sizes of problems for 
the computational study. The five sizes are (l2 x: 12), 

(24 x; 24), (48 x 48), (72 x 72), and (96 x 96), To explain 
the significance of the problem size, for example, let us 
consider the problem of (24 x 24) size. These numbers indicate 
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that the planning horizon for scheduling is 24 time periods 
and there are 24 different types of shifts available for assign- 
ing the workers. The motivation behind considering these- 
particular problem sizes is as follows. Normally, most of the 
service organizations such as telephone exchanges, hospitals-, 
police and fire departments work roim,d the clock. Depending 
on the demand, they can have a different shift starting every 
two hours, one hour, half an hour or quarter hour in a day* 

Thus, these particular sizes are considered. 

As we have already discussed in Chcpter I, a shift is 
differentiated from the other shifts by its start time or by 
its length or by the number of breaks or by the timing of 
breaks or by the combination of some or all of these. For a 
given problem, we differentiate two shifts by their start time 
in the problems considered. The cycle length, number of blacks 
of consecutive ones, q, and the timing of breaks are similar 
for all the n shifts. 

Within each size of the problem, we considered the shifts 
with different values of q. The values of q considered are 
2, 4, 6, and 8. For a given q and cycle length, we again consi- 
dered two cases. In the first case, the cycle is symmetric 
about its middle elements and in the other case the cycle is 
non— symmetric . Altogether, 31 problems, 16 symmetric and 
15 non-symmetric problems are considered* Tables 4,1 and 4,2 
respectively provide the details of the symmetric and non- 
symmotric problems considered. 
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A cycle is an (m x l) column vector. Given a cycle, we 
can generate the A matrix by taking the (m x l) column vector 
and shifting it cyclically one position dovaiward to obtain 
successive columns. Thus, the A matrix has the form 

^1 

^ 2 . 

« 

and is generated by the column vector (a^, a^^). 

Tt Is possible that not all the n shifts may be present 
in practice. As the number shifts becomes smaller, the problem 
becomes easier to solve. For all the problems are considered, 
all the n shifts are assumed to be present. 

The demands during different time intervals were generated 
randomly by taking the range as 100 to 150. 

In reality, there may be some pattern in the demand. 

To see how the performance of our algorithms is effected when 
there IS certain pattern in demand, we have considered the 
following five demand patterns and solved two problems of 
(A6 X 48) size and (72 x 72) size with a = 8. The demand 
patterns considered are as shown in Fig, 4,1, 




LEVEL 



SINUSOIDAL 


r" 


STEP 



4 J Various typ«s of demand 
patterns. 
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Here the demand is constant during all the time periods. 

Here, the demand increases linearly from start time 
period to the end time period, 

^ this case, the demand remains constant durissg first 
half of the time periods and then shoots up suddenly and 
remains at that value during rest of the time-periods* 

^j-P-P'SP.idals In xhis case the demand pattern is like a sine wave. 

Here, the demand pattern is like a sine wave during 
each half of the total time periods. 

In Tables 4,1 and 4,2, we have given certain number for 
each problem. In the subsequent tables and discussion, these 
problems are referred to by their numbers. 

Wo can observe that the CSS problem is completely speci- 
fied, If we are given the problem size, cycle and the range of 
demands . 

Computational Performances 

t -m mu <iki « M mmtdU m m in nn n m Mr t. s- jc.-. 

We, now give the computational perfoimiance of our algo- 
rithms and Bartholdi ^s algorithm in terms of accuracy and 
CPU times. 

In Table 4.3, we give the oboective function value obtained 
and the CPU times taken by Bartholdi ^s method, LP-based heuristic 
algorithm and the branch and bound algorithm for the symmetric 
problems. The total number of nodes enumerated in branch and 
bound algorithm is also given in the table. 
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It is evident from the table that the branch and bound 
algorithm was able to obtain optimum solution to all the 
problems. The maximum number of nodes enumerated are 53 for 
the problem 14, This problem has also taken the maximum 
amount of CPU time of 274.4 seconds. 

The objective function values obtained using LP-based 
heuristic are observed to be very close to the optimum value 
in all the problems. The deviation from the optimum solution 
Is less than 0.5 percent. Tho CPU time taken is much lesser 
when fomparod to branch and bound algorithm. This algorithm 
is observed to perform far better than Bartholdi "s method in 
terms of accuracy especially vdien q is more. 

It was observed that for none of the problems, the LP- 
based heuristic solved more than 2 LP'^s. In addition, the size 
of the second LP solved v;as often less than 10 percent of the 
original problem size. 

Bartholdi ^s method was observed to be often taking as 
many additional units as given by his bound. It is better only 
in terms of the CPU time taken. 

Similar data for the non-symmetric problems is given in 
Table 4.4. 

The performance of branch and bound is somevdiat discourag- 
6 of the 15 problems were solved optimally by 


ing in this case. 
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enumerating less than 25 nodes. In those cases where the 
algorithm did not terminate after enumerating 100 nodes are 
marked with a star. For these problems we usede - approximation 
wibh e = 1, In this case, we terminate the problem if the 
incumbent solutions [z] satisfies the following conditioni 

'2' is less than or equal to (Z + e ). 

The LP-based heuristic and Bartholdi'" s method behaved in 
the same way as tney did for the symmetric problems. 

Problems 8 and 12 were solved for each of the five demand 
patterns described in the earlier section and the results are 
given in Table 4.6. In these cases, the LP-based heuristic 
algorithm performance is same as for the earlier problems. The 
branch and bound algorithm was able to solve problem 8 optimally 
for all the demand patterns. Problem 12 could not be solved 
when the demand pattern is constant or sinusoidal. When 1 unit 
of e - approximation is applied, these tvro cases were also 
solved. From the nodes enumerated, we can observe that the 
problem becomes more difficult to solve when there is certain 
pattern in demand. 

The performance of the dual fractional cutting plane 
algorithm was quite erratic. The problems upto the size 
(72 X 72) were tried on this algorithm. Out of the 23 problems 
tried only 7 problems were solved and the others did not ter- 
minate after appending 150 cuts. The results of those problems 
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which were solved are given in Table 4.7. The algorithm could 
not solve problem sizes more than (48 x 48). In those cases, 
it could solve the problem, the number of cuts generated are 
very loss in number. 

4 , 4 Cpnclud^g^ 

In Chapter I, we observed that the CSS problem with 
intermittently available reso\jrces is more difficult to solve, 

¥o gave the statement of the problem and reviev/ed the litera- 
ture available for the CSS problem. We have developed three 
algorithms to solve the CSS problem more accurately. 

The LP-based heuristic algorithm developed in Chapter II 
makes uso of an improved version of the rounding-off procedure 
suggested by Bartholdi et al [ 5 ] and solves a sequence of LP'^s 
to got a good feasible solution to the CSS problem. 

Vfo developed a branch and bound algorithm in Qiapter III. 
The salient features of this algorithm areS (l) It makes use 
of the LP-based heuristic algorithm at every node to get a good 
feasible solution, ( 2 ) Sensitivity analysis was used to signifi- 
cantly reduce the storage and computational time requirements, 
and ( 3 ) the criterion used for source row selection makes prun- 
ing more effective. We have also developed a dual fractional 
cuttin/' plane algorithm in Chapter III. We proposed two types 
of cuts, one cut is based on objective function value and the 
other cut is Gomory'^s fractional cut. 
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Finally, we have done extensive computational study on 
three algorithms and Bartholdi* s algorithm. We 
observed that the performance of the LP~hased heuristic algo- 
rithm was quite encouraging. The deviation of the solution 
obtained using this algorithm from the optimum solution was 
not more than 2 units generally, ¥e observed that not more 
than 2 LP s were needed to be solved for any problem in the LP— 
heuristic algorithm and the size of the second LP solved is often 
less than 10 percent of the original problem size. The perfor- 
mance of the heuristic was invariant with the pattern of demand. 

The performance of fche branch and bound algorithm is quite 
encouraging. Symmetric problems of the size (96 x 96) with q = 8 
were solved optimally in a few minutes of CPU time. We have 
observed that problem becomes more difficult to solve when the 
cycle is non-symmetric or there is certain pattern in demand. 

The performance of cutting plane algorithm was observed 
to be erratic. 

We have considered the CSS problem when all c^ s are 1, 
Further scope for research in this area would be to extend the 
ideas developed in this thesis to consider the case when c^ s 


are different. 



Table 4.1s Symmetric Problems, 
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Table 4.3° Computational results for symmetric problems. 


U1 

tu sn 

H 

o 

G n • 

'i-i G 
O "H 
S 

• cd cd 
o x 
S J ffl 


K>HC\1C\I CM <i- U-\ <j- 



a 

O CO 

rH 

• 

cr\ 

• 

rr\ 

♦ 

cr\ 

* 

in 

« 

• 

H 

» 

CD 

• 

o 

• 

0 

rn 

• 

r—l 

• 

m 

« 

-d- 

N 

IS 


.cd 

w e 

o 

o 

o 



i>- 

CM 

M3 

ITS 

<1- 

in 

CM 

d- 

d" 

H 

in 

a 


rJ P- 





r 1 

(H 

CM 


H 

VD 

m 

H 

d* 

n- 

CO 

d- 

“ 

<■< -p^ 












H 


CM 


CM 

• 

l 

.to 


















.o 

•H 

o. 

iH 



in 

o- 

LA 

CTN 

m 

<1- 

LA 

VjO 

m 

cn 

d' 

CM 

CO 



.a- 

• 

• 

• 

• 

• 

• 

p 

• 

• 

• 

• 

• 

• 

* 

m 

* 

m 

J:=J 

O 4^, 

r 

O 

o 

o 

o 

m 

Lr\ 

rn 


m 

cn 

cn 

cn 

A 

IN- 

(Ts 

o 


)1 d 4-J 
'■H 

, W 

,o , 

<S .H 
i-H cd 
jH 

, jO -d 

P ifd o 
A.?P P 
P ,i^ -P 
- '&Si 


HH HHK^^-dLO 


ldCM<t CMrHPH©<l rdK^l>-IAIS-<r 

oo OK^L^^K^<l‘HIr^o^^^^Lr\'MOco^- 


5 fQ -P 
i,CD I 'H 

Id sd d 


cr\HO a^H^o^o^H[NCO(^JlNlr^POO 
jCTiOO l^~O^O^I^-O^O^O^OH^dCMK^^d 


iH io 
iP H 
{lO -”P 

K 


3 *H o 
^ 
r 

iw \ Cl) a> 

>ip ll 

p 

:o2 H 

p 

0 

, i o d 
t iJrd O 

1 ;p p 

» dp 

f ^ Si 

! .p S 


mcMoo mino CM CO (T^<t-roMDo hcm 
CTvO o CO a^o^®c^o^c^^o HrohorAro 


<}• CTirmiriCM-dHO t>-Pcomt>-LncM 
oocri(iiOcr\HooHrAK^K^-cfiA 
■ct-d- 


-ct-d-dprapro CMCMCMCMPCOVOP 
CM CM CM<f‘-d<J-<t' t^C^OtACTvCliCriCA 

<t■<^“<^COOOCOCX) C\JC\iCMC\lV£>UDU:)VjD 

CM CM cM^f^h'-d’^i" c^{>-{N-{>-cr\cr\0'>cn 


HCMrOi^j* (J\ O mUD 


H H H H H H H 



Table 4.4i Computational results for non-symmetric problems 


60 


m 

B ^ 
o 

'd • 
O -H 'd 


o ^ 
s w m 


i i 


pq I 

H 


jig o cgs 
3 bjOGrf 
I rH 
PQ -Pi 


; w , 

^ o , 3 
j 0) j ca 
;C0 . H O 
H ! P ^ I 
S .P o-pij 

J H 1, CU -H 0) 

;j ^P'-i -P S 

; a j I 

i H 4 cn , 

IiEH « i| 

j£ |h I 

ii '-p ^ ii 

,iq p ’ 


tl 


J CD O) 

m s 


HHK^OOOOmOOO 
CMCXJOOOHOOO 
H H H H H H 


/\ /\ /\ 


/\ A /\ 


o o o 
o o o 

H H H 
AAA 


<j- <}- 


^ LniAHHLnCMCr\OJCOlACOOOl>- 

oo<rir\i>-K^coLocriOcricr\co<i-H 

Pi — I O I — I -j" O 'd' ^ C\1 O I — I 

r’HH <J-K>,CI VOCOCTi 

r^/N/\ /\/\/\ /N/XA 


<r<i- LPiLPiiAvo'P'PLnHuq o-moouD 

• • «•««••••• •••« 

OOOK^<t-<ft'^^C^)C^^0^<^00<J-OO 

rH 1 — I P CM C\J <!■ CD <j" 


<f<r(^C7\OK^CDCMO--d-Ln[>-C7CCJ>VD 

OOOCMf^^-J-CMO^IpOOPOPOOc 

1 — I I — 1 I — 1 CM I'O <j" LD^ ro 




I :« I ! 

« 0) IIP P-i 

Ip !§ S g I 

i> 

j-H ; I 

ijP !• W 

J C » M 

ipip 


P-ci-VDOOOCDtO|POOPCOCMiM 

oooDoa<j\Chcr\coocr.oo<rpcMP 


o ] 

^ - 

H 


VO 

H 

O 

O 

CO 

•4 

H 

CM 

iH 

H 

00 

A 

4 

-P i! 

CX) 

CO 

0) 

Ch 


(T\ 

00 

O 

a\ 

O 

O 

4 

H 

CM 

1 — [ 

o u 

lA 


A 

A 

A 

A 

A 

4 

A 

4 

4 

4 

4 

4 

4 


!> ' w 

H 1 

Sp ,-h 

ii o 3 p 

i CD l,rH 
* CpS OP 
||P o 

10 Xi 

J i fnP 

J ! «J O 

11 iPq a 


p 

NS 

P 

CD S 


<1- CMOK^CO<t-VDCOCA~d-0<l-COCO CM 
CDCTiOCTiCrvOOOCrcPCM-cJ-CMK^K^ 


<1- <J- <t- CD CO CO CO 
CM CM CM -d" <t 


cmcmcmcmcocdcdco 

iPINC^C^CT'CTcOCCrv 


?SrSPSo?S?Sot^PSPN?S?SorSrS 


<t <1* <t CX) 00 CO CO 

C\] CM CM <t *4* <1* 


CMCMCM CMMOVOMDMO 

t^[>i>.i>^c^cr\cr\cr\ 


o • 

u o 

PM 12; 


OCOCTiOHCMK^^j-LAMOOCOCD^O H 
HHHCMCMCMCMCMCMCMCMCMCMIAIA 



I 


! U 


ftP! , 
2*2 i 

IN 


00 

cr\ 

H 

H 

lA 

A 

•4- 


•4' 

VO 

i ^ § 


cd "H I 

♦ 

• 

• 

• 

* 

• 

* 

# 

• 

• 

• 

• 

^ 1 CO 

40> 

. J 

! ^ 

CO S 



•4* 

VO 

ro 

0 

00 

H 

cn 

A 

A 

A 

j L 

.i OJ 

H 




(N 

rH 

OJ 

A 


C\J 

4- 

•4“ 

0 

A 


. PQ -P 


A 

H 

H 

A 

C\J 


(M 

CO 

C\j 

00 

00 

IN 


Sh 5 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

0 

P< 0 i 

A 

(N 

A 

CO 

cn 

CT\ 

0 

A 

Ch 

CO 

4 “ 

H 

g ^ 

cd 3 

H 

rH 

0 

rH 

- 4 - 

0 

* 4 * 

VO 

C>J 

A 

0 

H 

I 3 p 

■ - 


rH 

H 

H 


-4 

A 

OJ 


VO 

00 

cn 


! 

•H 

X 

o 

ft s:l 
ft o 
^ d -H 
-P , -P 
•H 8 cd 
^ to e 


I — ( rH 


'hS 


sipq 


H I 

X 1 

o 

■li ^ ^ ' 

d ftd , 

rt flj -H , 

-P , -P < 

•H J cd ij 

^ w s 


0 

0 

0 

0 

CTi 

0 

0 

0 

H 

0 

0 

0 

OH 

0 

0 

0 

rH 

0 

0 

0 


0 

0 

0 


H 

H 

rH 


H 

H 

rH 


H 

rH 

rH 


/\ 

/\ 

/\ 


/\ 

/\ 

/\ 



/\ 

/\ 


-H 

H X 

cd 1 o 

>|m P _ 

ll ft Sh 

^ 2*2 
os cd ’H 

•H i cd 40 -P 
-P Ij "H I cd 
o pq ^ w s 


rH 

0 

0 

VO 

4“ 

H 

0 

0 

H 

00 

OH 

CM 

A 

(J\ 

(T\ 

00 

0 

A 

0 

0 

4* 

rH 

C\J 

H 

A 

A 

A 

A 

4- 

A 

4“ 

4 - 

4- 

4- 

4- 

4" 


11 -h 1 

+> I of 

O PQ -P ^ J ^ 

oj 5 3 ft PJ J O 

nf Td o ft o j CTi 

^ 2 ^ cd •r 4 s f<^ 

o 3 40 ,40 t 
I "H I cd I 
!pq ^ oj g ! 


0 

^0 

*vO 

A 


*0 

*0 

H 

*00 

*CM 

*CM 

CT\ 

cn 

00 

0 

A 

0 

0 

4* 

H 

CM 

H 

A 

A 

A 

4- 

A 

4* 

4" 

4* 

4* 

4* 

4- 


CO 

CO 

00 

00 

CM 

CM 

CvJ 

CM 

VO 

VO 

VO 

VO 

4* 

4- 

4- 

4- 

IN 

IN 

tN 

N 

(y\ 

a\ 

A 

C7^ 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

00 

CO 

00 

00 

CM 

CM 

CM 

CM 

VO 

VO 

VO 

VO 

4- 

4- 

4“ 

4“ 

N 

[N 

IN 

N 

A 

a\ 

A 



H c\i -4- m 0- 

CM CM CM CM CM CM CM 


CO cr\ o H I, 
CM CM to ! 



Table 4.6: Performance of heuristic and branch and bound algorithms 
for various demand patterns. 
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